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PHYSICAL BASIS FOR A CONSTANT LAG TIME 
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ABSTRACT 

We show that the constant time lag prescription for tidal dissipation follows directly from the 
equations of motion of a tidally-forced viscous fluid body, given some basic assumptions. They are 
(i) dissipation results from a viscous force that is proportional to the velocity of the tidal flow (ii) 
tidal forcing and dissipation are weak and non-resonant (iii) the equilibrium structure of the forced 
body is spherically-symmetric. The lag time is an intrinsic property of the tidally-forced body and is 
independent of the orbital configuration. 



1. INTRODUCTION 

The origin of tidal dissipation in gaseous planets and 
stars is, for the most part, an un-solved problem. Due to 
the extreme weakness of dissipation in these nearly per- 
fect fluids, even identifying the correct theoretical frame- 
work to address this issue has proven to be difficult. De- 
spite this, there has been a recent revival in tidal theory 
that has produced many promising results (Ogilvie & Lin 
2004; Arras 2004; Wu 2005; Ivanov & Papaloizou 2007; 
Goodman & Lackner 2009; Weinberg et al. 2012). 

In order to study the long-term evolution of orbits that 
are shaped by tidal dissipation, a parameterization of the 
strength of tidal dissipation is often performed. One of 
the most common parameterizations is to incorporate a 
tidal lag time, where the forced response lags behind the 
equilibrium tidal deformation by some fixed value. Hut 
(1981) utilized this prescription in order to derive the 
orbit-averaged evolutionary equations for a wide array 
of orbital configurations. His results are widely used in 
many different astrophysical systems (e.g., Eggleton et 
al. 1998; Wu & Goldreich 2002; Wu & Murray 2003; 
Fabrycky & Tremaine 2007; Willems et al. 2010; Wu & 
Lithwick 2011; Socrates et al. 2012) 

A benefit of utilizing the constant time lag parameteri- 
zation is its ease of implementation. That is, the constant 
time lag model may serve as a useful tool when inferring 
the strength of tidal dissipation of similar objects that re- 
side in vastly different orbital configurations, which may 
in constraining actual theories of tidal dissipation (see 
e.g., Socrates et al. 2012). 

In order to further this goal, we attempt to understand 
the underlying assumptions behind the constant time lag 
model of Hut (1981). In his analysis, Hut approximates 
that the induced quadrapolar tidal deformation as two 
point particles, where the line joining them spans the 
diameter of the forced body and lags behind the line 
joining their center of mass and the distant source re- 
sponsible for the tidal field. Furthermore, the mass of 
the point particles is chosen such that the amplitude of 
their resulting quadrupole moment is equal to that of 
the equilibrium tide of a perfect fluid body. A deter- 
mination of the correspondence - if any - between the 
picture outlined above, to the fluid dynamics of a tidally 
forced body is the subject of this work. 

In Sj2]we demonstrate that, from the equations of mo- 



tion, of a non-resonantly weakly forced, spherically sym- 
metric body, the tidal response lags the equilibrium value 
by some fixed amount of time, in the event that the dis- 
sipation is weak. Furthermore, we show in f}2] that this 
constant lag time is an intrinsic property of the tidally 
forced body, in that it only depends upon its internal 
structure. A brief discussion and summary are given in 



2. ASSUMPTIONS AND DERIVATION OF THE CONSTANT 
LAG TIME MODEL 

2.1. the tidal interaction 

Consider a spherical, self gravitating fluid object 
(star/planet) which is weakly perturbed by a slowly vary- 
ing external gravitational field of the form 

U T (r,9,<f>,t) = -J2r e Y em (0^)V em (t). (1) 

em 

If, for example, the perturbation arises from a point mass 
m pcr orbiting the object at distance D(t) with angular 
coordinates 9' (t) , 4>' (t) we have 

lm ^ ' 

To leading order in perturbation theory, the interaction 
potential is given by (cf. Newcomb 1962) 



Hi = / d 3 xp£-VU T 



(3) 



where J d 3 x is taken over the volume of the forced body, 
p(x) is its unperturbed fluid density and £(x, t) is the 
Lagrangian displacement field. 
We may write the interaction energy as 



Cm 



(4) 



Cm 



where qe m is the multipole moment (cf. Press & Teukol- 
sky 1977) 



(t) = J d 3 xp$-Vr e Y t 



im V"i ' 



(5) 



Note that for I = 2, the \l/f m 's have dimensions of fre- 
quency squared. 
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The rate of energy transfer E between the orbit and 
and the forced body is given by 



E 



f d 3 xp£-VU T = -J2 

J Ira 



(6) 



Secular orbital evolution is entirely determined by the 
relationship between the the tidal potential, represented 
by the ^m's, and the corresponding induced multipo- 
lar moments, qem's. That is, the gravitational potential 
induced by the perturbation of the body's mass distri- 
bution is fully determined by the g^ m 's. While the fluid 
response £ may be complicated, only its multipolar mo- 
ments affect secular orbital evolution. 

2.2. equation of motion, the equilibrium tide and 
hiqher-order corrections 

We neglect rotation when considering the equilibrium 
structure as well as the dynamics of the fluid perturba- 
tions. Furthermore, we ignore the effects of rotation in 
the limit of non-resonant forcingQ The tidal problem 
may be expressed as 

£ +C£+D£= -VU T = -J2 ( 7 ) 

where C is an Hermitian operator that is responsible for 
the restoring force (Lynden-Bell & Ostriker 1967) and 
D is a time-independent differential operator that leads 
to dissipation. In the co-ordinate system of the forced 
body, the spatial dependance of the tidal forcing is given 
by the vector 



E f m (x) = v/r £m i 



(8) 



which for £ — 2 has dimensions of length. 

If the time-dependance of the forcing is slow such that 
the inertia of the fluid is small, then 

S < C • £ (9) 

and in the weak friction approximation 

D • £ <C C • £. (10) 

Given this ordering, we may approximate the solution 

as 

£(x, t ) = ^)-C- 1 .[D.£ (0) +£ (0) ] (11) 
and the equilibrium tide solution £^ is given by 

£<°> (x, t) = -J2 *im(t) Cr 1 • E £m (x). (12) 

2.3. a single lag time for a given £ 

The interaction energy, given by eq. 21 and the en- 
ergy transfer rate, given by cq. [51 only depend on the 
multipole moments qe m . When computing the multipole 
moments <# m 's by inserting eqs. 1111 and 1 1 21 into eq. [5l 
coefficients of the following form are encountered 



<i>S f * m .M-S Pr 



(13) 



1 Throughout, 'non-resonant' refers to the fluid oscillations of 
the forced body. 



where M = C" 1 or M = C ^D C 1 . By assuming that 
C and D are rotationally invariant, which is equivalent 
to assuming the equilibrium structure of the forced body 
is spherically symmetric, then 



(14) 



It follows that we can express a given multipole moment 
to leading order as 

= ~ n - Pi 2 (15) 

where a single lag time T£ for a given £ is 



and is given by 



Pi 



III 



(16) 



(17) 



Though the integrals above contain the azimuthal quan- 
tum number m, spherical symmetry requires invariance 
under rotations and consequently, independence of m. 

The term oc qf^ does not lead energy or angular mo- 
mentum transfer since it contribution is a full derivative 
in time, which does not accumulate over long time-scales. 
The leading order terms of the multipole deformation 
qem(t) that are responsible for apsidal precession and sec- 
ular orbital evolution are thus given by 

Urn (t) = qfi (t) - n $1 (t) ~ qW (t-n) (18) 

where we assumed that Tg is small in comparison to the 
characteristic time in which the tidal potential varies. 
Sole consideration of the quadrupolar £ — 2 response of 
the expression above is equivalent to eqs. 2 & 3 of Hut 
(1981), which together serve as the starting point and 
underlying assumption of his analysis. Furthermore, by 
inserting the above expression for the qi m into eq. [H we 
may write the secular energy transfer rate as 



E : 



E 



nqi 



(0)« 



(19) 



Again, by restricting to I = 2, the relation above is equiv- 
alent to eq. 40 of Eggleton et al. (1998), which serves 
as the starting point of their derivation for the secular 
equations for orbital evolution. 

The lag time th is completely determined by the equi- 
librium structure of the forced body. That is, T£ is de- 
termined by p(x), C and D all of which, under the as- 
sumptions previously mentioned, only depend upon the 
equilibrium structure of the object in question. Similar 
conclusions can be deduced from the analysis of Willem 
et al. (2010). Those authors considered the case of a 
spherically symmetric radiative star where dissipation re- 
sults from thermal diffusion and turbulent viscosity. 

2.4. constant density equilibrium structure 

As an illustrative example, consider a constant density 
equilibrium structure. Reisenegger (1994 and references 
therein) points out that for such an idealized system, 
Si m is an eigenfunction of C, with eigenvalue tog. Fur- 
thermore, and by construction, E^ m is responsible for 
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the entire multipolar response from the tidal accelera- 
tion — VC/t oc ai rn . In this case, the equilibrium tide 
becomes 

AO) _ ^im{t) „ 

Sim ~2~ a irn- (20) 

The displacement field B^ m is fundamental mode of 
the forced body, which is sometimes referred to as the 
/—mode or the Kelvin mode. S^ m has no radial nodes 
and its period of oscillation is, essentially, the free-fall 
time at the surface. 
In this limit, the lag time ti becomes 

where j£ is the damping rate of the H£ m 's and 

Pe="f- (22) 

Each Si m need not be an eigenvector of D in order for eq. 
|2"T1 to be correct. As in the more general case, the only 
requirement on D for producing a single constant time 
lag ti is for it to be invariant under rotations. Finally, 
for a constant density equilibrium structure, the tidal 
problem is equivalent to a set of decoupled forced damped 
harmonic oscillators with equal constant lag times. 

3. DISCUSSION AND SUMMARY 

The constant time lag model of Hut (1981) is perhaps 
the most widely used prescription for parameterizing the 
strength of tidal dissipation. Hut states that the most 
attractive feature of the constant r model for tidal dis- 
sipation is its simplicity. Perhaps this is true. However, 
as we have shown, the constant r model for tidal dis- 



sipation follows from some very basic physical assump- 
tions. Namely, the tidal forcing non-resonant, the forced 
body is a spherically symmetric fluid and the dissipation 
as well as the tidal forcing, is weak. In fact, these as- 
sumptions are well- approximated in many astrophysical 
environments and can be, in principle, tested with ob- 
servations of stars (cf. Dong et al. 2012) and extra-solar 
planets. 

We demonstrated that the lag time is an intrinsic 
property of the tidally-forccd object. In other words, two 
identical objects placed in vastly different orbital config- 
urations possess identical values for ti as long as the 
underlying assumptions of a constant Tg remain valid. 

For example, Socrates et al. (2012) show that in order 
for a Jupiter analogue to undergo high-e migration, the 
required lag time must be at least ten times stronger than 
that inferred from the Jupiter-Io interaction. Therefore, 
either the assumptions that led to a constant lag time are 
either incorrect or, for example, another element must be 
added to the theory of high— e migration. 
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